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ABSTRACT 


A rational methodology has been developed whereby three-dimensional 
sources and sinks may be placed along tne major axis of a class of ovary 
ellipsoids so as to minimize normal velocity and to calculate as exactly 
as possible the tangential velocity, pressure distribution, and the 
body shape. For this purpose the strength and position of the singular- 
ities and the position and number of the control points were optimized 
through the use of the method of least squares and the Automated Design 
Synthesis optimization technique. The results have shown that the 
previous methods are far from satisfactory and the use of two types of 
optimization in the determination of the strength and position of the 
Singularities yields the desired body shape and flow characteristics 
with excellent accuracy. A comprehensive computer code has been 
developed to enable one to calculate most of the practically 


Significant body shapes. 





TABLE OF CONTENTS 


CE ONC et == 2-22 eeee === ------------------------_-_-_.--.- 13 
ea smsg Diese == === -=------ -- ---------------5=-----.-..-- 16 
Semen Shen ICAL DESCRIPTION ------~----------~--~-~--------------_- 20 
ee nefenNAlnySiS AND RESULTS -------------------------------- 24 
A. USE OF EQUAL NUMBER OF SINGULARITIES 
iD CORAL, its sai eee ae 24 
Ee a SKOOEOE WEAST SQUARES =------------------------.-.--- 26 
ev ROPmaM I ZAunON MEMHODOLOGY -----------------------~----- 29 
SCT @RG Md) meee == ------------=--+--~--=-+--=--+---------- 35 
Se et Oy Weise = === ---------------- == --------.--------- 36 
Bee ee eERENCES ---~--~~~--~~~--~----------------- Sears eke oF 
tee CGO CUTER PROGRAM =--------------------~----~--=--------- 117 
SeeEE MECC EUMER PROGRAM =--------------~-~-~--~---~--~--~------=- 140 
INITIAL DISTRIBUTION LIST ------------------------------------------- 160 





Sul, 


Se). 
on. 
6a. 
6b. 
GG. 


6d. 


be. 
Sire 
7a. 
D. 
Ge, 


mis Or fb tGURES 


Ovary Ellipsoid -------------------------------------------------- 38 
Velocity Components ---------------------------------------------- 38 
Theoretical Tangential Velocity for Quarter Length of Body ------- 39 
Theoretical Tangential Velocity for Half Length of Body ---------- 40 
Placement of Singularities and Control Points -------------------- 4] 
Normal Velocity for the Case of 4/4 ------------------------------ 42 
Tangential Velocity for the Case of 4/4 -------------------------- 43 
Theoretical and Calculated Tangential Velocities 

for the Case of 4/4 -------------------------- -------------------- a+ 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Case of 4/4 ------------------------ 45 
Pressure Coefficient for the Case of 4/4 ------------------------- 46 
Body Shape for the Case of 4/4 ----------------------------------- 47 
Normal Velocity for the Case of 10/10 ---------------------------- 48 
Tangential Velocity for the Case of 10/10 ------------------------ 49 
Theoretical and Calculated Tangential Velocities 

for the Case of 10/10 ------------------------- 2-2-2222 een ne ---- 50 
Difference between the Theoretical and Calculated 

Tangential Velocities for Case of 10/10 -------------------------- 5] 
Pressure Coefficient for the Case of 10/10 ----------------------- eZ 
Body Shape for the Case of 10/10 --------------------------------- 53 
Normal Velocity for the Case of 20/20 ---------------------------- 54 
Tangential Velocity for the Case of 20/20 ------------------------ Sys 


Theoretical and Calculated Tangential 
Velocities for the Case of 20/20 --------------------------------- 96 





7d. 


/e. 
abi 
Sai. 
8b. 
Sc. 


8d. 


ee 
On. 
oa. 
9b. 
He . 


oar 


Je. 
Of. 


10a. 
WB. 
Oc . 


10d. 


10e. 
Chas 
lla. 


Difference between the Theoretical and Calculated 


Tangential Velocities for the Case of 20/20 ---------------------- 57 
Pressure Coefficient for the Case of 20/20 ----------------------- 58 
Body Shape for the Case of 20/20 --------------------------------- 59 
Normal Velocity for the Case of 4/11 ----------------------------- 60 
Tangential Velocity for the Case of 4/11 ------------------------- 6 | 
Theoretical and Calculated Tangential Velocities 

for the Case of 4/11 --------------------------------------------- 62 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Case of 4/1] ----------------------- 63 
Pressure Coefficient for the Case of 4/11 ------------------------ 64 
Body Shape for the Case of 4/11 ---------------------------------- 65 
Normal Velocity for the Case of 10/21 ---------------------------- 66 
Tangential Velocity for the Case of 10/21 ------------------------ 67 
Theoretical and Calculated Tangential Velocities 

for the Case of 10/21 -------------------------------------------- 68 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Case of 10/21 ---------------------- 69 
Pressure Coefficient for the Case of 10/21] ----------------------- 70 
Body Shape for the Case of 10/21 --------------------------------- TE) 
Normal Velocity for the Case of 20/41 ---------------------------- "eC 
Tangential Velocity for the Case of 20/41 ------------------------ He 
Theoretical and Calculated Tangential Velocities 

for the Case of 20/41 -------------------------------------------- 74 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Case of 20/41 ---------------------- 79 
Pressure Coefficient for the Case of 20/4] ----------------------- 76 
Body Shape for the Case of 20/41 --------------------------------- ve 
Normal Velocity for the Modified Case of 20/41 ------------------- 78 





Halabi. 
llc. 


| om 


Ile. 
it. 
ad. 
Zoe 
ZC. 


12. 


l12e. 
fez 
ieee 
ol. 
Ma. 
14b. 
14c. 


14d. 


14e. 
14f. 
fa - 
fob. 
[Dc . 


Tangential Velocity for the Modified Case of 20/41 -------------- 79 


Theoretical and Calculated Tangential Velocities 


for the Modified Case of 20/4] ---------------------------------- 80 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Modified Case of 20/41 ------------ 8] 
Pressure Coefficient for the Modified Case of 20/41 ------------- 82 
Body Shape for the Modified Case of 20/41 ----------------------- 83 
Normal Velocity for the Case of 20/40 --------------------------- 84 
Tangential Velocity for the Case of 20/40 ----------------------- 85 
Theoretical and Calculated Tangential Velocities 

for the Case of 20/40 ------------------------------------------- 86 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Case of 20/40 --------------------- 87 
Pressure Coefficient for the Case of 20/40 ---------------------- 88 
Body Shape for the Case of 20/40 -------------------------------- 89 
Initial Positions of Singularities and Control Points ----------- 90 
Final Positions of Singularities and Control Points ------------- oC 
Normal Velocity for the Case of 20/22 --------------------------- 9] 
Tangential Velocity for the Case of 20/22 ----------------------- a2 
Theoretical and Calculated Tangential Velocities 

for the Case of 20/22 ------------------------------------------- g5 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Case of 20/22 --------------------- 94 
Pressure Coefficient for the Case of 20/22 ---------------------- ole 
Body Shape for the Case of 20/22 -------------------------------- 96 
Normal Velocity for the Case of 10/12 --------------------------- oy 
Tangential Velocity for the Case of 10/12 ----------------------- 98 


Theoretical and Calculated Tangential Velocities 
for the Case of 10/12 ------------------------------------------- 99 





a 


OF . 


Se. 


19¢c. 


ed. 


19e. 


Lame 


Difference between the Theoretical and Caiculated 

Tangential Velocities for the Case of 10/12 --------------------- 100 
Pressure Coefficient for the Case of 10/12 ---------------------- 101 
Body Shape for the Case of 10/12 -------------------------------- 102 
Normal Velocity for the Case of 8/10 ---------------------------- is 
Tangential Velocity for the Case of 8/10 ------------------------ 104 
Theoretical and Calculated Tangential Velocities 

for the Case of 8/10 -------------------------------------------- 105 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Case of 8/10 -------~--------------- 106 
Pressure Coefficient for the Case of 8/10 ----------------------- 107 
Body Shape for the Case of 8/10 --------------------------------- 108 
Tapered Axisymmetric Body Used in Analysis ---------------------- 109 
Theoretical Tangential Velocity for the 

Tapered Axisymmetric Body --------------------------------------- 110 
Calculated Normal Velocity for the 

Tapered Axisymmetric Body --------------------- 9 enone isle 
Calculated Tangential Velocity for the 

Tapered Axisymmetric Body --------------------------------------- 2 
Theoretical and Calculated Tangential Velocities 

for the Tapered Axisymmetric Body ------------------------------- ele: 
Difference between the Theoretical and Calculated 

Tangential Velocities for the Tapered Axisymmetric Body --------- 114 
Calculated Pressure Coefficient for the 

Tapered Axisymmetric Body --------------------------------------- Ike 
Calculated Body Shape for the Tapered Axisymmetric Body --------- 116 





TABLE OF ABBREVIATIONS AND SYMBOLS 


One half of the length of the major axis of an ellipsoid 
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IT. INTRODUCTION 


The development of reliable methods for the calculation of three- 
dimensional viscous flows on shapes of hydrodynamic interest requires 
the analysis of (a) the inviscid flow, (b) three-dimensional laminar 
and turbulent boundary layers, (c) flow separation, (d) post-separation 
flow and wake, and (e) the interaction between the viscous and inviscid 
flow regions. 

A critical survey of the present state of the art indicates that 
many aspects of the problem remain tentative or unresolved even though 
considerable progress has been made in recent years in improving the 
calculation methods for inviscid flow and boundary layers about three- 
dimensional bodies. Further progress in each of these areas is hampered 
largely by the lack of pertinent data from carefully conceived and 
executed experiments. To be sure, such experiments are difficult and 
very expensive. The availability of high-speed computers, on the other 
hand, has led to the development of novel calculation procedures which 
claim a level of generality that has far surpassed the reliability of 
the underlying assumptions and techniques. In the present study, 
attention is focused on the inviscid flow about axisymmetric bodies 
and its calculation through the use of exact and approximate methods 
of optimization to the desired degree of accuracy. 

The potential flow about an arbitrary body of revolution was first 
treated by von Karman [1]. He determined the potential flow around 


bodies of revolution at zero angle of attack by superposing a uniform 
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stream on a system of Tine sources distributed along the axis of the body. 
An equal number of coordinates (control points) were chosen on the body 
‘and the strengths of the sources were determined so that the zero stream- 
line passed through the chosen coordinates. 

Since the pioneering work of von Karman, very little original work 
has appeared in the literature. While most fluid dynamicists are familiar 
with the method of representation of a body with distributed sources (sinks) 
and surface panels, they generally do not recognize the limitations of the 
methods and provide no guidance as to how one can systematically approach 
the ideal solution. Often, broad statements and vague suggestions are 
made regarding the slenderness of the body, the number of the singularities 
(i.e., sources and sinks), and the magnitude of the differences between 
the calculated and exact solutions. 

Another approach to solving the inviscid flow problem for bodies of 
revolution is the surface singularity or surface panel method. This 
approach is equivalent to the solution of an integral equation (Fredholm 
integral of the second kind). However, the computational effort in the 
axial-singularity-distribution metnod of von Karman is a fraction of that 
of the surface singularity methods. Furthermore, the surface singularity 
methods are not necessarily more accurate than the methods of discrete 
or continuous singularities along the axis. 

It is, therefore, necessary to define more precisely the number of 
Singularities, their strengths and positions, the number of control points, 
the differences between the calculated and given body shape, the deviation 
of the normal velocity from zero on the body surface, the discontinuities 


or artificial oscillations in the tangential velocity and in the pressure 
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distribution on the body surface, and as to how one can assess the degree 
of accuracy of the calculation procedures and improve systematically the 
accuracy of the quantities calculated. This investigation deals with 
these questions and significantly improves the calculation of the incom- 
pressible potential flow about bodies of revolution at zero angie of 


attack. 





II. PREVIOUS STUDIES 


Von Karman [1] in 1927 used a method involving line sources clustered 
in the bow and line sinks clustered in the stern of an airship. The 
method of solution separated the bow computations from the stern computa- 
tions since the influence of one on the other was considered smal]. In 
this method a good comparison of airship pressures with experiment was 
obtained. However, for shapes more complex than airships one would have 
to consider sources and sinks along the whole axis of the body and one 
would not isolate discrete parts of the body for separate computations. 
In doing so, however, one discovers that the method is unstable or may 
fail unexpectedly. 

Hess [2] and Oberkampf and Watson [3] have critically examined the 
generalizations of von Karman's method. They have shown that the method 
produces a system of linear equations which is, in general, i]]1-conditioned, 
and requires very high computational accuracy in the construction of the 
coefficient matrix and in solving the equations. Oberkampf and Watson 
concluded that the method does not always produce reliable solutions for 
the flow around a specified body and that the conditions that the body 
Should meet in order to be represented by a system of axial l]ine sources 
are not clear. Karamcheti [4] states that the body should be slender 
and should not have any discontinuities in the slope of the meridian 
line. Numerical experiments by Oberkampf and Watson showed that such 
conditions are not sufficient. For example, the method gave a slight 


rippling in the velocity distribution as well as in the meridian 
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streamline of very slender Rankine ovals. Oberkampf and Watson attributed 
the rippling in velocity to the increased local effects of each source 
element. Ancther result of their study is that the method is sensitive 

not only to the shape of the body contour but also to the number of elements 
used to generate the body. QOberkampf and Watson attempted to calculate the 
potential flow about a sphere using an odd number of sources vice an even 
number. The results for an odd number of sources were very peculiar. 

While the zero streamline did pass thrcugh all of the specified coordinates 
of the sphere, the body generated had "holes" in its surface. That is, 
between the specified coordinates the zero streamline plunged into the nega- 
tive strength line sources on the axis and then reappeared so that it passed 
through the next specified coordinate. This behavior produced highly un- 
realistic normal and tangential velocities. It is hypothesized that this 
erratic behavior was caused by the fact that with an odd number of line 
sources one of the sources must overlap the center of the sphere to preserve 
symmetry and, consequently, prohibits any antisymmetric solution. Therefore, 
odd numbers of sources must be avoided for bodies that are symmetric about 

a plane normal to their major axis at the midpoint of their length. 

It is evident from the foregoing that von Karman's method does not 
always produce reliable solutions for the potential flow around a speci- 
fied body. The alternative is to use surface singularities or panels 

(2, 4]. The panel methods could calculate all the required viscous and 
inviscid flow properties without further interaction on the part of the 
user once the geometry of the body and the flow conditions are specified. 
However, such methods require far more computation time than any other 


inviscid flow method. - Accordingly, many investigators have concentrated 
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on representation of the inviscid flow through the use of discrete 
singularities (three-dimensional point sources and sinks) placed on 
the axis of the body. Once the surface velocity and pressure are cal- 
culated, the potential flow and the boundary layer are “patched together" 
to represent the real flow. It is important to note that in all of 
these calculations, regardless of the method of representation of the 
inviscid flow, it is assumed that the flow is essentially unseparated and 
the effects of viscosity are appreciable only in a very thin layer adja- 
cent to the body surface and in a thin wake downstream of the body. 

The ultimate goal of the use of discrete singularities along the 
axis of the body is to determine a flow which is, in some sense, a good 
approximation to the exact flow about the given body. If this approach 
is adopted, it is essential to know, however, the degree of approximation. 
This can be checked in a number of ways after one sclves for the strength 
and/or the position of the singularities. One laborious way of checking 
the accuracy would be to increase the number of singularities, solve 
again for the strengths, positions, and velocities, and compare the two 
solutions for convergence. This process has three essential drawbacks. 
First, it can lead to inefficient use of the computer facilities. The 
computer time required for the solution increases rapidly with the number 
of singularities. Thus, the cost for solution depends very strongly on 
the number of the singularities used. If the number of singularities 
is increased until little change in the local velocities on the body is 
observed, then certainly the final and by far the most expensive computa- 
tion (i.e., the last one) is redundant. Second, as discussed in the 


foregoing, such a process may not converge when distributed singularities 
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are used. Third, and perhaps most important, there is nothing in this 
process which yields an insight into how one might improve the strength 
sjaeposition of the singularities. In other words, it is difficult to 
determine if a more efficient computation would result if one region of 
the axis had more singularities and control points per unit length. 

An alternative method of checking the accuracy of the solution is to 
evaluate the norma! velocities on the exact surface of the body (i.e., 
on the shape specified rather than on the one represented by the zero 
streamline) at points other than the control points. If the body is 
nonporous, as assumed here, then these normal velocities should be 
identically zero. It is easily shown that the magnitude of the normal 
velocities on the exact body is, in fact, a proper measure of the in- 
accuracies in the whole flow field. Thus, the remnant normal velocity 
distribution at the end of a particular step may be used to improve the 
calculations and hence the calculated flow field. 

In view of the simplicity of the use of discrete sources and sinks 
along the axis of the body as compared to the distributed axial singular- 
ities or surface singularity distributions, and the tremendous saving in 
the amount of memory storage and numerical calculation, it was considered 
justifiable to undertake an extensive study of the body representation 
by discrete singularities. In this study, improvement is realized by 
allowing the strength and position of the singularities and the number 
and positions of the control points, respectively, to achieve their 


Optimum values. 
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TIT. MATHEMATICAL DESCRIPTION 


In this section, the basic equations for a discrete axial distribu- 
tion of three dimensional sources and sinks combined with uniform flow 
are derived. In addition, the exact tangential velocity for an ovary 
ellipsoid is obtained for comparison with that obtained numerically. 

Stokes' stream function for an ambient flow of unit velocity past an 
axisymmetric body with N sources (sinks) of strength Q; along its axis 


is given by [5] 





where ro = at + (x - a) at the point (x,w) in the flow field. Evidently, 


w = 0 corresponds to the enclosed body (see Fig. 1). 


The velocity components u and v are given by 





N Nomar 
usa Seles m. <4 (2) 
ans is] r 
N 
5. Ls w 
es a ox u ae (3) 
T=] Yr 


where m. = Q;/4r. 

The body is assumed to be defined either by a function w = f(x} or 
by a discrete set of points. In either case, there is sufficient infor- 
mation to calculate the normal and tangential components of the velocity 


along the body. In fact, from Fig. 2 one has 
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u =v cos@ - u sin@ (4) 

V, =v sind + u cos@ (5) 
where @ is defined by 

tan@ = df/dx (6) 


It is evident from the foregoing that the accurate specification of the 
stream function determines the accuracy of the tangential velocity and 
the existence of non-zero normal velocities along a non-porous wal] 
indicates the error between the calculated and ideal solution. 

The numerical experiments have been carried out with ellipsoids of 
various a/b ratios (see Fig. 1). Thus, it was necessary to derive the 
exact expression for the tangential velocity so that the accuracy of 
the numerical method may be properly assessed. 


The transformation 


oti = C COSNt va 
yields 

Kee Coss | COST (8) 
and 

Cees ime. SIN} (9) 


Thus, & = EO denotes an ellipse, in the meridian plane, as 


2 Ze 


Xx W) 
a te Ci (10) 
ic coshé ) (c sinhe )“ 

whose semi-axes are 

a=c coshé . io = Siite sinhé | (1) 


Note that b/a = tanhé | determines EO 
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It is relatively easy to show that the stream functicn is given by 





: 2 e 

Ww =U cf sinh°e sinen pes (coshé + sinhcé > Ln tanh 5) sac (V2) 
2 2K 2 
in which 
2 
ea 8 + 2 i ae) 
Ce ahs 1) iia aaa 
The velocity components are given by [6] 
-__l oy -_!| sy 

We Jo n ee een) = Ge 14) 
where 

ce (15) 
and 

2, ‘ 1 

Jee ft (ec) f (cz) (16) 
and 

mee) = € Coshc and f (rc) =e ecuiiae ep) 
Equations (16) and (17) yield 

J = ae (Cosh 2a COs) (18) 
V2 

Thus, the velocity components on the body (i.e., for & = ee ) are 

/2 cosn bo So 
te = ee | SINE = aE, (ctnhé ouiiieesen tanh zee (19) 

Ycosh 2 .7coséan Ke : 


AE 





’ and 





2 2 
(1 + coshé -Ln tanh 5) (20) 
Y/cosh 26 7cosen Ke c 


It can be shown that Eq. (19) is identically zero since it represents 
the normal velocity on the body. The tangential velocity is calculated 
from Eq. (20) for representative values of a/b and is shown in Figs. 3a 
and 3b as a function of the normalized distance from the forward stagna- 


tron point. 
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IV. NUMERICAL ANALYSIS AND RESULTS 


A. USE OF EQUAL NUMBER OF SINGULARITIES AND CONTROL POINTS 

As noted in connection with the discussion of previous investigations, 
the use of the singularity methods requires the selection of the position 
of the singularities along the axis and of the control points on the 
body. A convenient but arbitrary selection of boundary points and 
Singularity locations in this procedure is to position the singularities 
directly below the boundary points in one-to-one correspondence [6]. 
Evidently, no special criteria 1s provided for the spacing or number of 
the singularities. In any case, this procedure results in a set of linear 
equations which can be solved through the use of standard matrix reduction 
techniques. 

The first example chosen to illustrate the technique and the problems 
associated with it is an ovary ellipsoid (also called a prolate spheroid, 
generated by rotation of an ellipse about its major axis). The ellipsoid 
had a slenderness ratio of a/b = 6.0. An even number of singularities 
(N.) was chosen and they were equally spaced along the major axis of 
length 2a. An equal number of control points No was placed directly 
above the singularities on the ellipse as shown in Fig. 4. 

Figures 5a through 5f show the results obtained with four singularities 
and four control points. Figure 5a shows the normal component of the 
velocity (V) calculated along the upper half of the ellipse. The con- 
dition of w = 0 is exactly satisfied at the control points. However, 


large non-zero normal velocities between the control points show clearly 
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the inadequacy of the body representation by four equally spaced singular- 
ities and control points. Figure 5b shows the tangential velocity (V,) 
along the body. Theoretically, one would expect a smoothly increasing 
tangential velocity profile. Thus, the large oscillations in vee are a 
further indication of the failure of only four singularities and control 
points, as presently positioned, to adequately represent the body. Figure 
5c shows a comparison of the theoretical and calculated V. as a function 
of the normalized distance from the forward stagnation point while Fig. 
5d depicts the difference between tneoretical and calculated tangential 
velocity along the first half cf the body. Apparent!y, the calculated 
Vy oscillates about the exact velocity profile. In Fig. 5e the pressure- 
coefficient C,, [C, = (p - p.)/(0.5 u°)], and in Fig. 5f the body shape 
resulting from the 4/4 case are presented. It is evident from the fore- 
going that an ellipsoid of a/b = 6.0 cannot be successfully represented 
by a small number of equally spaced singularities and control points. 
Thus, it is necessary to explore first the effect of increasing the 
number of singularities in the representation of the body and flow 
characteristics prior to embarking on a more ambitious investigation of 
the effect of using nonuniform singularity distributions. 

In anticipation of improvement of the flow characteristics with an 
increased number of singularities and control points the cases of 10/10 
and 20/20 were investigated. Figures 6a through 6f and Figs. 7a through 


7f show the results for the cases 10/10 and 20/20, respectively. 


Rte convenience, the designation N./N., e.g., 4/4, will be used 
hereafter to refer to the number of singularities and control points 
used in a particular example. 
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Evidently, the difference between the theoretical and calculated tangen- 
tial velocity, as well as, all other errors in the calculated flow char- 
acteristics have been reduced by increasing the number of singularities 
and control points. Nevertheless, the representation of the flow char- 
acteristics is far from satisfactory and one cannot use the resulting 
Stream function to predict the boundary layer characteristics over the 
body. It is, therefore, necessary to explore other methods which will 
minimize the number of the singularities as well as the error between 


the predicted and calculated flow characteristics. 


B. THE METHOD OF LEAST SQUARES 

The use of an equal number of singularities and control points resulted 
in a deterministic set of linear equations and in making w exactly equal to 
zero at each and every control point., However, one could use a larger 
number of control points tnan singularities. This will obviously result 
in an over-determined set of linear equations if one attempted to render 
w = 0 at all control points. Recognizing the impossibility of doing so, 
One can, instead, minimize w at all control points through the use of a 
Suitable minimization technique. Thus, one can make the error in w nearly 
zero at a larger number of control points in lieu of making it exactly 
zero at a fewer number of control points. 

Let us consider N. Singularities and Ne control points. Tne stream 
function at the control point j due to the contributions of all singular- 


ities is given by (see Eq. (1)), 
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If w were not exactly zero at the point j, then the error in i would be, 


Ng X . X - 


R. = m. aces) _ a y 4 (22) 
1=] r —_ 


Then the sum of the square of the errors at all control points becomes, 


Gee 
j=l 


Z 


N 
—- 2 


(23) 


in which x is a function of the strengths of Ne number of singularities 
only. The total error may be minimized througn the use of Gauss' method 
[7] by writing the partial derivative of x with respect to each m, 


eamal to zero, 1.é., 





fa (26) 
r 
[b] = [5 04] (27) 
and 
m = m. (28) 


In the foregoing, i varies from | to Ne and j, from | to Ne 
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The least squares method described above was applied to the ellipsoid 
with a/b = 6.0 for the cases 4/11, 10/21, and 20/41. As before the 
Singularities were equally spaced along the major axis. The control 
points were likewise located on the body with equal spacing. The results 
are shown in Figs. 8a through 8f for the case of 4/11, in Figs. 9a 
through 9f for the case of 10/21, and in Figs. 10a through i0f for the 
case of 20/41. Tne comparison of these figures with those cited earlier 
(j.e., with those having identical letter designations such as 5a and 8a, 
etc.) shows that the increase in the number of control points does in 
fact improve the calculated tangential velocity and reduce normal velocity. 
One can also observe some improvement in the body representation as 
evidenced by a comparison of Figs. 5f and 8f or of Figs. 7f and l10f. 

It is also clear that an ellipsoid of slenderness ratio of a/b = 6.0 
Cannot be adequately represented with four singularities even with a 
larger number of control points. Furthermore, even with a significantly 
larger number of singularities and control points (as in the case of 
20/41) the method of least squares does not yield a satisfactory repre- 
sentation of the flow characteristics. This leads one to the conclusion 
that optimizing the strengths of the singularities alone is not sufficient 
to adequately model the flow about the body. 

The realization of this fact leads to the hypothesis that the locations 
of the singularities and the control points may be less than optimum. To 
test this hypothesis, one could judiciously position both the singularities 
and the control points in the areas of the body where the error is large 
(j.e., large normal velocities) and then compare the calculated flow 


Characteristics with those obtained in the case of evenly spaced 
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Singularities and control points. Previous calculations have shown that 
the error in all flow characteristics is largest where the radius of 
curvature of the body is relatively small (i.e., regions near the fore 
and aft stagnation points). With this in mind, the two singularities 
closest to the center of the body in the case of 20/41 were moved to 

the new positions of x/b = + 5.71429 while the remainder of the singu- 
larities were kept in their original evenly spaced positions (see Fig. 1). 
The four control points closest to the center of the body were moved to 
new positions directly above x/b = + 5.90476 and x/b = + 5.80952 while 

the remainder of the control points were kept in their original positions. 
Figures lla through 11f show the results obtained in this manner. A 
comparison of these Figures with Figs. !0a through l10f shows that Ue 

has been significantly reduced and there is greater agreement between 

the theoretical and calculated tangential velocities as a result of the 
repositioning of the singularities and the control points. Thus, one 
must devise a rational method with which the singularities and the control 
points can move to their optimum positions as the singularities continue 
to adjust their strengths, as in the previous cases, through the use of 


the method of least squares. 


C. DUAL OPTIMIZATION METHODOLOGY 

In the following, a method is described whereby the position and 
strength of the singularities and the number and position of the control 
points are progressively optimized so as to minimize the error in the 
prediction of the flow characteristics about a given body. For this 


purpose the following steps have been developed: 
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1. Decide on the appropriate number of singularities, Nes for the 
body of given slenderness ratio (two singularities per unit length are 
recommended on the basis of the experience gained with the least squares 
method). 

2. Position the singularities with equal spacing along the major 
axis of the body. 

3. Select twice as many control points, Nee as Singularities and 
place them on the body contour with equal horizontal spacing. 

4. Determine the strength of the singularities through the use of 
the method of least squares, as described previously. 

5. Calculate the sum of the absolute values of the normal velocities 
between singularities N. and Ne ay for a arom 1 te No as well as, in the 
regions between the forward stagnation point and the first singularity 
and between the last singularity and the rear stagnation point. This 
yields Ne) sums corresponding to the (N_+2) intervals. In addition, 
calculate the total sum of the absolute values of the normal velocities 
along the full length of the body (hereafter referred to as = }). 

6. Initially, place one control point per interval on the body 
contour. This requires (N. + 2) control points. The remaining 
LN. - (N 


using the relative magnitude of the sum of the absolute values of the 


+ 2)] control points are apportioned among the intervals 
normal velocities in each interval as a weighing factor. Specifically, 
the extra control points are assigned in accordance with 

x 


N [= =e 


Aer : acy (23) 
ony |G) 
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where ee represents the number cf control points at the end of a given 
iteration and k varies from 1 to (N42). As progressive iterations of 
this process reduce the normal velocities along the body, less than the 
original number of control points may be necessary to adequately represent 
the body and the flow characteristics. Equation (29) allows for this and 
will reduce the total number of control points as needed. However, the 
number of control points will never be less than (N.+2). 

7. Place in each interval the apportioned number of control points 
with equal horizontal spacing. 

8. Determine the improved positions for the Me Singularities through 
the use of the Automated Design Synthesis optimization program [8]. 

9. Repeat steps 4 through 8 until no further reduction is realized 
mi, 

The complete computer program based on the steps described above is 
presented in Appendix A. 

Several representative calculations have been performed to demonstrate 
the effectiveness of the dual optimization methodology. The first case 
concerns the ellipsoid with the slenderness ratio of a/b = 6.0. Twenty 
Singularities and 40 control points were chosen in accordance with the 
Suggestions made in Steps # 1 and 3. At the end of the first iteration 
(j.e., Steps #1 through 9) the total number of the control points needed 
was reduced to 36. Figures 12a through 12f show the results in graphical 
form. Clearly, the normal velocities are fairly large and the difference 
between the theoretical and calculated tangential velocity is not yet 
satisfactory. The computer program carried out a total number of 28 


iterations at the conclusion of which X was reduced to its minimum value. 
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After the final iteration, the total number of the control pcints was 
reduced to 22. The control points and singularities became unevenly 
spaced and moved towards the fore and aft regions of the body. Figures 
13a and 13b show, respectively, the initial and final positions of the 
Singularities and control points. Figures 14a through 14f show that the 
difference between the theoretical and calculated tangential velocity is 
reduced to about 0.10xU very near the stagnation points and to less than 
about 0.04xU along the remainder of the body. Furthermore, the calculated 
velocities over the central half of the body are almost identical with 
those predicted theoretically. 

It was noted earlier that the dual optimization program should begin 
with two singularities per unit length. To validate the significance of 
this suggestion and to provide another test case for the computer code, 
the case of 10/20 was considered for the ellipsoid of slenderness ratio 
of a/b = 6.0. The dual optimization program performed as expected and 
at the end of 22 iterations achieved the minimum 2 . In the course of 
the optimization, the number of the control points was reduced to 12. 
The final results are shown in Figs. 15a through 15f. A comparison 
of these figures with Figs. 14a through 14f shows that the use of less 
than two singularities per unit length does not produce flow character- 
istics as good as those calculated in the case of 20/22. Thus, the 
use of about two singularities per unit length is considered optimum. 
Obviously, it is always possible to use a larger number of singularities. 
However, the additional computing expense is not commensurate with the 


minimal improvement achieved in the flow characteristics. 
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To evaluate the ability of the computer code to deal with ellipsoids 
of other slenderness ratics, a test case with a/b = 2.0 was run. 
Initially, 8 singularities and 16 control points were equally spaced 
along the body. Only three iterations were required to achieve the 
minimum value of x. The results, presented in Fig. 16a through l6f, 
show that the difference between the calculated and theoretical velocities 
has been reduced to about 0.04xU along the entire body. Note also that 
the number control points required was reduced to 10 during the dual 
optimization process. 

To further establish the generality of the dual optimization method- 
ology, axisymmetric bodies without fore and aft symmetry have been con- 
sidered. For bodies of this type there is no simple boundary function 
in the form of w = f(x). The body shape is in general provided by the 
naval architect in accordance with the needs of the user. In the present 
Study such a body has been generated through the use of one source and 
four sinks. The strengths and the positions of the sources and sinks 
were such that the resulting body did not have fore and aft symmetry. 


The Stokes stream function is given by 


| 2s (X-x,) 
a (30) 
2 1=] 2 Zaye 
L(x-x,) + w,] 
where 

m, = 0.24 Nay 0.48891 

My = -Q.06 Xo = sya eletens) 

Mz = -0.06 X= 7.48891 

M4 = -0.06 Xq = 8.9889] 

Me = -0..06 Xp =10.4889] 


$8) 





The resulting body shape is Shown in Fig. 17. The exact tangential 
velocity calculated through the use of Eq. (30) is shown in Fig. 18. 
It will be compared later with that obtained numerically. 

The optimization process followed the steps outlined previously. 
Thus, in accordance with the recommendation that abcut two singularities 
per unit length be chosen, 20 singularities and 41 control pcints were 
selected and distributed appropriately along the body. [In only two 
iterations the number of control points reduced to 32 and £ acquired its 
minimum value. 

Figures 19a through 19f show the results after the final iteration. 
Clearly, the calculated and the exact body shapes are nearly identical 
and the normal velocity has been reduced to almost zero along the full 
length of the body. Also, the calculated and theoretical tangential 
velocities compare exceedingly well. [It should be nected that the un- 
dulations in the tangential velocity at the aft end of the body are due 
to the geometry of the body (as defined by Eq. (30)) and are not attri- 
butable to any theoretical or numerical instability. In fact, the 
occurrence of such undulations in the exact tangential velocity profile 
provided greater challenge for the dual optimization process. The 
computer program for the case under consideration is presented in 
Appendix B. 

It is evident from the foregoing that the dual optimization method- 
ology can be used with great confidence in the prediction of the flow 
characteristics about axisymmetric bodies with or without fore and aft 


Symmetry. 
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Ve) SCCHCEUSIONS 


The investigation described herein warranted the following conclusions: 

1. The existing methods for the representation of axisymmetric bodies 
and the flow about them require excessively large computer time and a 
great deal of foresight for the selection of the number and position of 
the singularities and control points so as to achieve satisfactory results. 
Furthermore, they provide no insight as to how the errors (e.g., normal 
velocity) may be progressively reduced. 

2. Through the use of the method of least squares and the Automated 
Design Synthesis optimization, together with a sufficient number of dis- 
crete singularities and control points, one can represent an axisymmetric 
body and the flow about it with excellent accuracy. The bodies are not 
required to have fore and aft symmetry. 

3. Extensive calculations with various types of axisymmetric bodies 
have shown that the use of two singularities per unit length is quite 
adequate to model the flow about the body. 


4. Numerous examples have been given and the results have been 


compared with those obtained theoretically. 





VI. RECOMMENDATIONS 


The following recommendations are made for the purpose of increasing 
the power of prediction of the dual optimization methodology. 

]. The existing computer program should be converted to an inter- 
active mode so as to simplify its use. 

2. the code should be enhanced to enable the user to minimize not 
only the strength and location of the singularities but also their 
number. For example, such a procedure could be incorporated into the 
code by deleting from the calculations the singularities whose strengths 
fall below a prescribed percentage of the sum of the absolute values of 
the strengths of all singularities. 

3. The improved code should be tested with other axisymmetric body 


Shapes. 
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Figure 1. Ovary Ellipsoid 





Figure 2. Velocity Components 
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Figure 4. Placement of Singularities and Control Points 
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